Parametric stiffness analysis of the Orthoglide 
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q \ This paper presents a parametric stiffness analysis of the Orthoglide. A compli- 

ant modeling and a symbolic expression of the stiffness matrix are conducted. This 
allows a simple systematic analysis of the influence of the geometric design param- 
eters and to quickly identify the critical link parameters. Our symbolic model is 
used to display the stiffest areas of the workspace for a specific machining task. Our 



Abstract 



o 

\& • approach can be applied to any parallel manipulator for which stiffness is a critical 

issue. 
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Usually, parallel manipulators are claimed to offer good stiffness and accuracy prop- 
erties, as well as good dynamic performances. This makes them attractive for in- 
novative machine-tool structures for high speed machining [1,2,3]. When a parallel 
manipulator is intended to become a Parallel Kinematic Machine (PKM), stiffness 
becomes a very important issue in its design [4,5,6]. This paper presents a parametric 
stiffness analysis of the Orthoglide, a 3-axis translational PKM prototype developed 
at IRCCyN [7]. 

Finite Element Methods (FEM) are mandatory to carry out the final design of a 
PKM [8]. However, a comprehensive three-dimensional FEM analysis may prove 
difficult, since one must repeatedly re-mesh the PKM structure to determine stiff- 
ness performances in the whole workspace, which is time consuming. Simpler and 
faster methods are needed at a pre-design stage. One of the first efficient stiffness 
analysis methods for parallel mechanisms was based on a kinetostatic modeling [9] . 
According to this approach, the stiffness of parallel mechanisms is mapped onto 



their workspace by taking into account the compliance of the actuated joints only. 
It is used and complemented in [10] to show the influence of the compliance of the 
prismatic joints as well as the torsional compliance of the links on the stiffness of the 
3-UPU mechanism assembled for translation [11]. It is shown that the compliance 
of the links reduces the kinetostatic performances in a large part of the workspace, 
compared to the stiffness model based on rigid links. Furthermore, the mobile plat- 
form can undergo small rotational motions because of the links' compliance, which 
departs from the expected translational kinematic behavior. 

The analysis presented in [9] is not appropriate for PKM whose legs, unlike hexapods, 
are subject to bending [12]. This problem is solved in [13], where a stiffness estima- 
tion of a tripod-based overconstrained PKM is proposed. According to this approach, 
the PKM structure is decomposed into two substructures, one for the mechanism 
and another for the frame. One stiffness model is derived for each substructure. The 
superposition principle allows one to join the two models in order to derive the stiff- 
ness model of the whole structure. The influence of the geometrical parameters on 
the stiffness is also briefly studied. An interesting aspect of this method is that it 
can deal with overconstrained structures. However this stiffness model is not gen- 
eral enough. A more general model was proposed in [14]. The method is based on a 
flexible-link lumped parameter model that replaces the compliance of the links by 
localized virtual joints and rigid links. The latter approach differs from that pre- 
sented in [13] on two main points, namely: (i) the modeling of the link compliances 
and (ii) the more general nature of the equations allowing the computation of the 
stiffness model. 

In this paper, the method proposed in [14] is applied to the Orthoglide for a paramet- 
ric stiffness analysis. A symbolic expression of the stiffness matrix is obtained which 
allows a global analysis of the influence of the Orthoglide's critical design parameters. 
No numerical computations are conducted until graphical results are generated. This 
paper is organized as follows: first the Orthoglide is presented. Then, the compliant 
model is introduced and the stiffness model is computed. Analytical expressions of 
the components of the stiffness matrix are obtained at the isotropic configuration, 
clearly showing the influence of each geometrical parameter. Finally, given a specific 
simulated machining task, it is shown how the general stiffness expressions allow one 
to easily display the stiffest subvolume of the Orthoglide's workspace. 



2 Compliant modeling of the Orthoglide 

2. 1 Kinematic architecture of the Orthoglide 

The Orthoglide is a translational 3-axis PKM prototype designed for machining 
applications. The mobile platform is connected to three orthogonal linear drives 
through three identical RP a R serial chains (Fig.l). Here, R stands for a revolute 
joint and P a for a parallelogram-based joint. The Orthoglide moves in the Cartesian 




Figure 1. The Orthoglide (a) Kinematic architecture and (b) Prototype 

workspace while maintaining a fixed orientation. The Orthoglide was optimized for 
a prescribed workspace with prescribed kinetostatic performances [15]. Its kinematic 
analysis, design and optimization are fully described in [15]. 

2.2 Parameters for compliant modeling 

The parameters used for the compliant modeling of the Orthoglide are presented on 
Fig. 2 and in Tab. 1. They correspond to a "beam-like" modeling of the Orthoglide 
legs' links. The foot has been designed to prevent each parallelogram from colliding 
with the corresponding linear motion guide. Three revolute joints are added, one 
on each leg (see Fig. 2), because the stiffness method used does not work with an 
overconstrained Orthoglide. This does not change the kinematics. 




Figure 2. Geometric parameters of the leg Figure 3. General model for a flexible link 

(a) Flexible beam (b) Virtual rigid beam 



2.3 Compliant modeling with flexible links 

In the lumped model described in [14] , the leg links are considered as flexible beams 
and are replaced by rigid beams mounted on revolute joints plus torsional springs 
located at the joints (Fig. 3). Deriving the relationship between the force F and the 
deformation y(x), the local torsional stiffness k can be computed: 



EIy"{x) = F(L - x) 

Ely(L) = FL 3 /3 
y(L)/L = FL 2 /3EI 
k = FL/9 
-»• k = 3EI/L 

If the Orthoglide leg actuator is locked, then one leg can withstand one force F and 
one torque T (Fig. 4), which are transmitted along the parallelogram bars and the 
foot. For a compliant modeling that uses virtual joints, it is important to understand 
how external forces are transmitted, and what their effect on the leg links is. Eight 




Figure 4. Forces transmitted in a leg 



virtual joints are modeled along the Orthoglide leg. They are described in Tab. 2. 
The determination of all the virtual joint stiffnesses is not detailed here for brevity. 
However, they are derived based on the same principles used to calculate the torsional 
stiffness above. 



Parameter 


Description 


Values 


Lf 


Foot length, see Fig. 2 


150mm 


hf 


Foot section sides 


26mm 


bf 


Foot section sides 


16mm 


1 h 12 


Foot section moment of inertia 1 




h f .b) 
L h ~ 12 


Foot section moment of inertia 2 




I f0 = h f .b f (hj + b})/12 


Foot section polar quadratic moment 




A 


Angle between foot axis and actuated joint axis, 
see Fig. 2 


45° 


d 


Distance between parallelogram bars, see Fig. 2 


80mm 


L B 


Parallelogram bar length, see Fig. 2 


310mm 


S b 


Parallelogram bar cross-section area 


144m?7i 2 


P 


Rotation angle of the parallelogram 




e 


See Fig.2 




Table 1 



Geometric parameters of the Orthoglide and dimensions of the protptype 



The actuated joint is assumed to be much stiffer than the virtual joints. The leg links 
compliances modeled in Tab. 2 were selected beforehand as the most significant ones. 
Indeed, selecting only the most significant compliances plays an important role in 
reducing the computing time required to derive the stiffness matrix symbolically 
(Par. 3). The kinematic joints' compliances are not taken into account because our 
purpose is to determine the links compliance influence only. Angle (3 is a parameter 
that depends on the Cartesian coordinates. 



3 Symbolic derivation of the stiffness matrix 



In this section, the derivation of the Orthoglide stiffness matrix — based on the 
virtual joints described in the previous section — is conducted with a stiffness model 
that was fully described in [14]. Therefore, the description of the model will only be 
summarized here. Fig. 5 represents the lumped model of a leg with flexible links. The 
Jacobian matrix Jj of the ith leg of the Orthoglide is obtained from the Denavit- 
Hartenberg parameters of the ith leg with flexible links. This matrix maps all leg 
joint rates (including the virtual joints) into the generalized velocity of the platform, 
i.e., 

■T 



J;0j = t where 0, 
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Virtual joints i 
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Virtual joints i 



Figure 



k\ — kact 
translational 
stiffness of the 
prismatic 
actuator 



EI 



- Lf 

Foot section 
rotation due to 
torque T 




3EI 



Foot bending 
due to force F 




Eb 

Parallelogram 
bars tension/ 
compression due 
to force F 




2EI 



k 3 

Foot bending 
due to torque T 




k 



ES B d 2 cos(fl) 
10 - 2li 

Differential 
tension of 
parallelogram bars 
due to torque T 




Gl 



So 



-4" Lf 

Foot torsion 
due to torque T 




Table 2 
Virtual joints modeling 



Figure 5. Flexible leg 

is the vector containing the 11 actuated, passive and virtual joint rates of leg i and 
t is the twist of the platform. The P a joint parameterization imposes 8i 7 = —8i 7bis , 
which makes 9 ir and 9 irbis dependent. 9 ir is chosen to model the circular translational 
motion, and finally Jj is written as 
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in which e^. is the unit vector along joint j of leg i and is the vector connecting 
joint j of leg i to the platform reference point. Therefore the Jacobian matrix of the 
Orthoglide can be written as: 



Ji 
J 2 
J 3 



One then has: 

30 = Rt with R = [I 6 I 6 I 6 ] T and t 



6 being the vector of the 33 joint rates, that is 6 = [0 1 6 2 6 3 ] T . I 6 stands for the 
6x6 identity matrix. Unactuated joints are then eliminated by writing the geometric 
conditions that constrain the two independent closed-loop kinematic chains of the 
Orthoglide kinematic structure: 



Ji#i = J 2 #2 and J\6i = Js0 3 



(2) 



From (2), one can obtain AO' = BO" (see [14] for details), where is the vector of 
joint rates without passive joints and 6 is the vector of joint rates with only passive 
joints. Hence: 

e" = b'ao' 

Then a matrix V is obtained (see [14] for details) such that: 



e = ve 



(3) 



From (1) and (3) one can obtain: 
JV6>' = Rt 



(4) 



As matrix R represents a system of 18 compatible linear equations in 6 unknowns, 
one can use the least-square solution to obtain an exact solution from (4): 

t = (R T R) 1 R T JV0 

Now let J' be represented as J' = (R T R) X R T JV. Then one has: 

t = J'O' (5) 



According to the principle of virtual work, one has: 

r T e' = w T t 



(6) 



where r is the vector of forces and torques applied at each actuated or virtual joint 
and w is the external wrench applied at the end effector, point P. Gravitational 
forces are neglected. By substituting (5) in (6), one can obtain: 



J /T w 



(7) 



The forces and displacements of each actuated or virtual joint can be related by 
Hooke's law, that is for the whole structure one has: 



K,A0' 



(8) 



with Kj : 



A 
A 
A 



and A = diag ( k act , ^, 2fSa ^»f~(fl 



AO' only includes the actuated and virtual joints, that is by equating (7) with (8): 

KjA0' = J ,T w 

Hence AO' = Kj 1 J /T w. Pre-multiplying both sides by J' one obtains: 

J' AO' = J'Kj^w (9) 



Substituting (5) into (9), one obtains: 

d = J'Kj^w 

with d = tAt. Finally the compliance matrix n is obtained as follows: 
In the Orthoglide case we obtain: 
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«n 





K 24 K 25 










«11 


«34 «35 


«36 


«i4 


K 24 




K44 K45 


K 46 


«15 


^25 




K45 K55 




V^16 


«26 


«36 


^46 K 56 


«66 / 



(10) 



And the Cartesian stiffness matrix is: 



1 t'Tn- 1 



K = k~ = (J'Kj- 1 J' 1 ) 



4 Parametric stiffness analysis at the isotropic configuration 



In this section, we study the influence of the geometric parameters on the stiffness 
of the Orthoglide at the isotropic configuration, since this configuration provides a 
good evaluation of the overall performances [15]. Another interest is that the stiffness 
matrix is then diagonal which makes it easier to analyze. 



4-1 Simple symbolic expressions 



At the isotropic configuration, n is diagonal and the symbolic expressions of the 
components are simple. This is convenient because it is then possible to invert 



k within a Maple worksheet and then analyze the symbolic expressions of the com- 
ponents of matrix K. We have: 

K = diag (tf a , K a , K a , K b , K b , K b 

where K a is the torsional stiffness and K b is the translational stiffness. 

E 1 

Ka = 2L B 2L p (78fe 2 +cos 2 A(45fe 2 -33fe^ Kb = [ . L B , 4L^. sin 2 A ( U ' 

Sgd 2 ^ 5hfb 3 f (b 2 f +h 2 f ) kact + 2S B E + Ehjbf 



Analyzing the Orthoglide's stiffness at the isotropic configuration allows us to ma- 
nipulate simple and meaningful symbolic expressions that are easy to interpret: this 
is the purpose of the following subsections. 



4-2 Qualitative analysis of K a and K b 



The expression — lh f b s {b 2 +h 2 ) ~ shows the influence of virtual joints 3, 4 



By inspection of the symbolic expression of K a a few observations can be made: 

• Young's modulus E appears at the numerator, which makes its influence easy to 
understand: when E increases, K a increases, which is in accordance with intuition; 

• The term J^Jj- shows the influence of virtual joint 10 (differential tension of par- 
allelogram bars). When the bar length Lb increases or when Sb decreases, K a 
decreases which is also in accordance with intuition. K a decreases when d in- 
creases, which is a less intuitive resuliPl: 

2L p (786 2 +cos 2 A(45/i 2 -33f> 2 )) 
b){b 2 +h 2 ) ~ 

and 5 (foot bending and torsion). K a decreases when Lf increases, which is not 
surprising. The degrees of hf and bf in the numerator and denominator of K a tend 
to prove that the rotational stiffness increases with hf or bf, which is in accordance 
with intuition. The influence of A depends on the sign of (45/i? — 336j). 

Similarly, by inspection of the symbolic expression of K b one notes: 

• The term j^— shows the influence of the prismatic actuator: it is not surprising that 
the translational stiffness increases when k act increases. The term 2 s b e sri ows the 
influence of virtual joint 8 (parallelogram bars tension/compression): K b increases 
when Sb or E increase, and decreases when L B increases, which is in accordance 
with intuition; 

• The term BhFb} shows the influence of the foot related virtual joints (ten- 
sion/compression and bending): when A increases, with A e [0 rr/2] 0, then sin A 



1 Note that should d increase above a certain limit, other links compliances previously 
ruled out as less significant may then need to be taken into account. 

2 If A > 7r/2 the foot does not anymore "move away" the parallelogram from the prismatic 



increases and consequently Kb decreases. According to intuition, increasing Lf 
decreases Kb, while increasing hf or bf increases Kb- 



4-3 Quantitative analysis of K a and Kb 



As we have seen, the qualitative analysis of K a and Kb provides interesting infor- 
mation on the influence of the geometrical parameters on the rotational and trans- 
lational stiffnesses. Quantitative information about the parameters' influence on the 
Orthoglide's stiffness can also be obtained from the symbolic expressions by study- 
ing the consequences of a - 100/+200% variation of the parameters on K a and K b . A 
variation of -100% corresponds to a zero parameter, while +200% corresponds to an 
extreme increase. Such a wide range of variation gives a global picture of the param- 
eter's influence. The initial values of the parameters used for the computation are 
given in Tab. 1 and correspond to the dimensions of the prototype of the Orthoglide 
developed at IRCCyN. Parameters k act and E are considered constant because our 
analysis is restricted to geometrical parameters only. We choose E = 7.10 4 Nmm~ 2 
(aluminum) and k ac t = 10 5 Nmm" 1 . The stiffness of the actuated prismatic joint 
depends on many parameters (mechanical components, electrical motor power, con- 
trol). The chosen value is a commonly used one, however it is still much stiffer than 
the virtual joints, which is in accordance with our assumptions. 

In order to clearly show the relative influence of each parameter, we are going to 
superimpose several curves on a same chart. Each curve represents a ratio R Ka ^ 

a initial 

(resp. K Kb ^ ), in which t is the percentage of variation of one of the parameters 

initial 

(Lf, bf, hf, X, L B , Sb or d), while the other parameters remain at their initial value, 
and K a . nitial (resp. K hmitial ) is the initial value of the torsional (resp. translational) 
stiffness when the parameters are at their initial value. Obviously, all R Ka ^ (resp. 

a initial 

K b (t) 



^initial 



) curves cross when t = 0%. 



For example let us replace each parameter in the symbolic expression of K a by its 
initial value except Lf. A one variable analytical expression K a (Lf) is then obtained: 

In this expression, let us replace Lf by Lf. nitial (1 + t). A new expression K a (t) is 
obtained: 

„ . . 0.56 x 10 9 

KJt) = ; r 

° W 150(1 + *) 

where t represents the percentage of variation of Lf. K a (t = 0) gives the value for 
K a . nit . al . We assume that t varies from —100% to +200% as explained above. All 



actuated joint, which is one of its main functions (i.e. avoiding collisions between the 
actuator and the parallelogram); furthermore we must have A > to avoid interference 
between the foot and the actuated prismatic joint. 



K a {t) / K ainitial curves obtained for all parameters are superimposed on a same chart 
so as to compare the parameters relative influence. 

• Quantitative analysis of K a 

Figure 6 shows the influence of the parameters on K a . Lb, d and Sb have little 
influence compared to Lf, hf, bf and A. 



K a (t)/K a . nitial : most influent parameters K a (t)/K a . nitial : least influent parameters 



K a (X) is a maximum (52% increase) when A increases by 100%, i.e. when A = 7r/2. 
This result can also be obtained through observation of the symbolic expression of 
K a : indeed, the initial values of hf and bf (hf = 26, bf = 16) make (45/t? — 33b"j) 
positive. Therefore, the denominator of K a will be a minimum when A = 7r/2. 
Moreover, when A = tt/2, the torque T that is transmitted by the leg no longer has 
a component along the axis of virtual joints 3 and 5 of the foot (Fig. 7). This is a 
physical explanation for K a (X) being maximum when A = ir/2. 



Figure 7. Only virtual joints 4 of the foot is affected by T when A = tt/2 

Furthermore, K a increases more with bf than with hf for a same variation. Conse- 
quently, for a given foot weight increase, the torsional stiffness benefits more from 
an increase of bf than from an increase of hf. From a designer's point of view, this is 
valuable information. If the foot length Lf increases, K a decreases since in this case 
the foot and torque related stiffnesses k 3 , A; 4 , k 6 decrease. Conversely if Lf decreases 
then K a increases tremendously. 




t variation (%) 



t variation (%) 



Figure 6. Influence of the parameters on K ( 



a 




Finally, we observe that when d, Sb, hf or bf tend towards zero, then so does K a . 
This can be deduced from the symbolic expression of K a , but also tends to a physical 
interpretation: if the foot or the parallelogram bars tend to have a very small cross 



section, or if the parallelogram tends not to be able to support any torque (when d 
tends towards zero), then the whole mechanism loses its torsional stiffness. Though 
hf and A play important roles in K a , the two most important parameters are Lf 
and bf. 

• Quantitative analysis of Kb 

Figure 8 shows the influence of the geometrical parameters on K^. We can observe 
that Lb and Sb have little influence compared to Lf, hf, bf and A. Kb(X) is a 
minimum (48% decrease) for a 100% increase of A, i.e. when A = tt/2. This conclusion 
can be reached through the observation of the symbolic expression of Kb. indeed, 
we can see that the denominator of Kb will be a maximum when A = n/2. 



Lb 

I 

'~Sa 

-50 50 100 150 200 

t variation (%) t variation (%) 

Kb{t) / Kb initial : most influent parameters Kbit) / it al : least influent parameters 
Figure 8. Influence of the parameters on Kb 




From the symbolic expression of Kb, one can also infer that if A decreases, then the 
denominator will decrease and consequently Kb will increase. This was the opposite 
case for K a . For a 100% decrease of A, Kb will be a maximum: 14.4 times its initial 
value. This has a physical interpretation: when A = 0, the virtual joint 2 is no longer 
affected by the force F that is transmitted by the leg (Fig. 9). 




Figure 9. Only virtual joint 3, 4 and 5 of the foot is affected by F when A = 



Kb is also a maximum (14.4 times its initial value) when Lf = 0. However, the 
physical interpretation is not the same. When Lf = 0, the stiffness of the virtual 
joints 2, 3, 4 and 5 tends toward +oo which makes them behave like infinitely stiff 
virtual joints, making the mechanism as a whole much more stiffer. One can also 
observe that Kb increases more with hf than with bf. This can be concluded from 



the symbolic expression of K b . Indeed, from Eq. 11, we have: 



Kh{kh hf) 0.00002537698413 + ^J 143 

Consequently, a 10% increase of hf will make the denominator of K b decrease faster 
than a 10% increase of bf. Finally, if S B , hf or bf tend towards zero, then K b also 
tends towards zero. This can be concluded from the symbolic expression of K b , but 
it also corresponds to the physical phenomenon that was explained for K a . 

The most important parameters for K b are A, Lf and hf. Parameters A and Lf have 
a similar influence: when they decrease, K b increases, and conversely. Parameter hf 
has the opposite influence: when hf increases, K b increases, and conversely. The 
symbolic expressions of K b as univariate functions of these three parameters are of 
great help at a pre-design stage to analyze the translational stiffness. 



4-4 Conclusions 



The analysis of the symbolic expressions of K a and K b at the isotropic configuration 
allows us to plot the most influent parameters in this configuration, and the way their 
variation influences the mechanism's stiffness. A global analysis must be conducted 
in the whole workspace to determine the global influence of the parameters. This was 
achieved in [17], with the determination of a line along which the stiffness analysis 
results hold for the whole workspace. Such a procedure is required to simplify the 
global stiffness analysis. However, as mentioned above in the case of the Orthoglide, 
analyzing the stiffness at the isotropic configuration can give a good overview of the 
performances. 

The use of simple symbolic expressions allows us to deduce helpful results in order 
to improve the Orthoglide's stiffness. However these modifications must be made 
while taking into account the technological constraints (collisions, interferences) of 
the prototype initial architecture. For example, if one sets A to zero in order to 
increase K b , then the offset between kinematic joints L 6 and L\ disappears. However, 
this offset aims at preventing the parallelogram from colliding with the prismatic 
actuated joint. Therefore it is not possible to set Lf or A to zero. It is better, either to 
only lower them and check how much the reachable workspace is then reduced, or to 
increase hf, or both. Conversely, if one wants to increase the collision-free workspace 
by increasing Lf while keeping K b constant, studying the simultaneous influence on 
the stiffness of Lf and hf or Lf and bf can then prove useful. One problem will be the 
foot weight increase that will require more powerful actuators to keep the dynamic 
performances at a similar level. We consider the issue of simultaneous variation of 
two parameters in the following section. 



5 Influence of the simultaneous variation of two parameters 



In this section, we study the influence on K a and K b of the simultaneous variation of 
two parameters Lf and hf, or Lf and bf, at the isotropic configuration. Analytical 
expressions of K a and K b as functions of two variables are deduced from Eq. 11. 
Figures 10 and 11 show plots of K a /K a . nitial and K b /K b . nitial when t hf (resp. t bf ) - 
which is the relative variation of hf (resp. bf) — and £5. — which is the relative 
variation of Lf — increase from to 100% or 200% when relevant. 




K a \thfi^Lf) /Ka initia i K a (t bf ,tL f ) /K a . niUal 

Figure 10. K a / K a . nitial as a function of hf, bf and Lf 



Figure 10 shows that increasing hf or bf allows us to compensate for the decrease 
of K a occurring when Lf increases. For example if Lf increases by 50%, hf must 
increase by 34% or bf must increase by 16% for K a to remain at its initial value 
K a . niUal . Regarding the dynamic performances (i.e. the foot weight increase), it will 
be more interesting to increase bf by 16%. 




K b (t hf ,t Lf )/ K b . n . Ual K b (t b f ,t Lj )f K bmiUal 

Figure 11. K b /K b . nitial as a function of hf, bf and Lf 

On Fig. 11, one can also observe that increasing hf or bf allows us to easily com- 
pensate for the decrease of K b occurring when Lf increases. If Lf increases by 50%, 
hf must increase by 48% or bf must increase by 245% for K b to remain at its initial 



value Kb initial . Therefore, it seems more judicious to increase hf rather than bf in 
order to compensate for the stiffness loss due to the increase of L/, because the foot 
weight increase is lower. This is a multi-criteria multi-parameters (Lf, hf, bf) op- 
timization problem: increasing the collision-free workspace while keeping the same 
stiffness, with a minimum foot weight increase. Integrating the symbolic expressions 
of the stiffness in multicriteria optimization loops could be an interesting extension 
of our work. 



6 Analysis of the tool displacements induced by external forces 



Another interesting use of the symbolic expressions of f% is to observe the tool 
compliant displacements when simulated cutting forces are applied on the tool. By 
multiplying these forces with the compliance matrix and analyzing the evolution of 
the compliant displacements obtained, as a function of the Cartesian coordinates, 
the stiffest zones of the mechanism's workspace can be determined. Thus, the global 
stiffness behavior is taken into account. As the simulated cutting forces correspond 
to a particular manufacturing operation, the stiffest zone will be specific to the 
application. The equations with which the stiffness matrix is computed are built 
using the principle of virtual work. Simulated cutting forces will then correspond to 
quasi-static conditions, which may not be realistic in some cases. In this section, a 
simple groove milling operation is simulated, which can be considered as a quasi- 
static operation. 

The symbolic derivation of the stiffness matrix K using the method described above 
was achieved with Maple software on a 1 GHz, 256MB RAM PC. The computation of 
K did not end within one day, which means that the components of matrix K, i.e. the 
Kij, are too large to be manipulated within a Maple worksheet. However, computing 
the components of the compliance matrix n, i.e. the k^, took 12 hours only. This 
resulted in symbolic expressions that remained relatively easy to manipulate within 
a Maple worksheet. Therefore we choose to analyze the Orthoglide's stiffness through 
the analysis of the symbolic expressions of the k^: the main idea is that when the 
increase, then the Orthoglide's stiffness decreases. 



6.1 Compliant displacements 



Vector w is the static wrench of the cutting forces applied on the tool during the 
groove milling operation along the y axis. We have: 



w = < 



T 
F 



> with F = [F x F y F Z ] T and T = [-F y h z F x h z 0} T (Fig. 12) 
The compliant displacements of the mobile platform are computed as follows: 



d = kw with d 




ft = \uj x ojy lu z ] t and V = [v x v y v z ] T 



The compliant displacements at the tool tip are then: 




V + ft x h z z 



6.2 Determination of the stiff est working zones for a given task 

With the symbolic expressions of the tool displacement, one can evaluate the tracking 
error along the groove path. Using the symbolic expression of the tracking error, a 
stiffness favorable working zone, i.e. a working zone in which the tracking error is low, 
can be determined. To simulate cutting forces during the groove milling operation, 
a High Speed Machining (HSM) simulation software is used [18]. Depending on the 
manufacturing conditions, this software provides the average cutting forces. The 
manufacturing conditions chosen for the groove milling are: 

• Spindle rate is N=20,000tr.min _1 ; 

• Feed rate V/=40m.min _1 ; 

• Cutting thickness is 5.10~ 3 mm; 

• The tool is a ball head of $ = 10 mm diameter with 2 steel blocks; 

• Manufactured material is a common steel alloy with chromium and molybdenum. 

The simulated cutting forces correspond to a HSM context, which is what PKM 
are F x = 215N, F y = — ION, F z = —25N. The above data allows us to simulate 
the tool compliant displacement along a groove path along the y axis (see Fig. 12). 
h z =100 mm corresponds to the tool mounted on the prototype of the Orthoglide. 
The tracking error is the projection of the tool compliant displacement in the plane 
that is perpendicular to the path. We specify one groove path with its coordinates 
(xt,Zt), and one point P with (x t ,yp,z t ) coordinates located along this trajectory. 



The paths are defined in a cube centered at the intersection of the prismatic joints, 
Xt and Zt vary within the interval [—73.65; 126.35]. We noticed that the maximum 
tracking errors were always located at one of the path ends. Figure 13 shows the 



The tracking error at point P is defined as Sp 





Figure 12. Component forces of groove milling operation 



tracking error along a groove defined with the coordinates (x t ,z t ) = (0,0). We 
can see that the maximum error occurs when y = —73.65, i.e. at one of the path 
ends. Depending on the coordinates (x t ,z t ), the maximal tracking error is located 
at y = —73.65 or at y = f26.35. Figure 6.2 shows the maximum tracking error for 
each groove path defined by its coordinates (x t , z t ). The results clearly show a zone 
in which the maximum tracking error is low. In this working zone, x varies within 
the interval [-73.65;0] and z varies within [50; 126.35]. It is difficult to find a physical 
explanation for this result. It depends on the cutting forces applied, their magnitude 
and direction, and on each virtual joint reaction to the wrench transmitted by the 
leg, which depends on the Cartesian coordinates. The information obtained, i.e., the 
lowest tracking error working zone, is, however, of great interest for the end-user in 
order to place manufacturing paths in the workspace achieving the lowest tracking 
error due to structural compliance. 




Figure 13. Tracking error along the groove Figure 14. Maximum tracking error along 
path defined by x = z = y axis groove paths 

Another use of the symbolic expressions of the compliant displacements would be 
the optimization of the geometric parameters to minimize the tracking error for 
specific cutting forces. This would mean optimizing a PKM design for a specific 
task. Our opinion is that it is better to look for global stiffness improvement as we 
did in the previous section. This way, optimization brings stiffness improvement to all 
potential manufacturing tasks. However, given a PKM design, it is very interesting 
to determine the stiffest working zones for specific tasks, as we did in this section. 



7 Comparison with a Finite Element Stiffness Model 

By comparing our stiffness model with a Finite Element Model (FEM) of the Or- 
thoglide prototype, we will now show that our rigid link model is reasonably realistic 
[17]. A FEM was implemented in LARAMA (LAboratoire de Recherches en Automa- 
tique et Mecanique Avancee, Clermont-Ferrand, France) as part of a collaboration 
within project ROBEA, a research program sponsored by CNRS (Centre National 



de la Recherche Scientifique). Due to space limitations, the modeling assumptions of 
the FEM are not detailed here. The FEM allows to calculate the variation range of 
diagonal elements Kt^i, Kt 2 ,2, Kt 3)3 of translational stiffness matrix Kt, based on 
a CAD model of the Orthoglide implemented in the finite element software ANSYS 
[19]. The results obtained are presented in Tab. 3 (deterministic approximations and 
variation ranges) at the isotropic configuration. Our objective is to compare these 
results to those obtained with our Rigid Link Compliant Model (RLCM). Stiffnesses 
are expressed in N.mm -1 . The numbers obtained from the FEM are comparable 





Kt M 


Kt 2 ,2 


Kt 3)3 


FEM 


RLCM 


FEM 


RLCM 


FEM 


RLCM 


Isotropic configuration 


3500 


2715 


3500/4000 


2715 


3500/4000 


2715 



Tabic 3 



Comparison of the RLCM and the FEM 

to those obtained from the RLCM. Even if deterministic values are not equal, this 
comparison shows that the RLCM of the Orthoglide is reliable enough for the pur- 
poses of pre-design. However, a more detailed FEM analysis and experimental results 
based on the Orthoglide prototype would be necessary to validate our RLCM. The 
main advantage of the RLCM is that it allows to spot critical links within the whole 
workspace much more easily and quickly than the FEM, because of the symbolic 
expressions of stiffness matrix elements. The RLCM is easier to use than a FEM at 
a pre-design stage. Once the RLCM is proved reliable enough, one can use it either 
to test alternative designs or choose manufacturing paths reducing the tool compli- 
ant displacement (tracking error or tracking rotational and translational compliant 
displacements) caused by structural compliance. Unfortunately, the FEM did not 
provide any results for the rotational compliance. This would be an interesting com- 
parison since it would allow a verification of whether or not the torsional stiffness 
of the mobile platform obtained with the RLCM is lower compared to that of the 
overconstrained Orthoglide prototype described in [15,16] and modeled in the FEM. 



8 Conclusions 

In this paper, a parametric stiffness analysis of a 3-axis PKM prototype, the Or- 
thoglide, was conducted. First, a compliant model of the Orthoglide was obtained, 
then a method for parallel manipulators stiffness analysis was applied, and the stiff- 
ness matrix elements were computed symbolically in the isotropic configuration. In 
this configuration, the influence of the geometric design parameters on the rotational 
and translational stiffnesses was studied through qualitative and quantitative anal- 
ysis. The analysis provided relevant and precise information for stiffness-oriented 
optimization of the Orthoglide. Then, the analysis of the simultaneous influence on 
the stiffness of two variable parameters was conducted. Such an analysis is very 
useful to take into account both stiffness and another performance criterion such 
as workspace volume or the maximal acceleration of the mobile platform. Finally, 



we used the symbolic expressions of the components of the compliance matrix to 
determine the stiffest working zone for a specific manufacturing task. The stiffest 
zone depends on the task and applied cutting forces. The parametric stiffness anal- 
ysis shows that simple symbolic expressions carefully built and interpreted provide 
much information on the stiffness features of parallel manipulators, which can be 
relevantly used for their design and optimization. 
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